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We discuss classical magnetotransport in a two-dimensional system with strong scatterers. Even 
in the limit of very low field, when lu c t <C 1 (uj c is the cyclotron frequency, r is the scattering 
time) such a system demonstrates strong negative magnetoresistance caused by non-Markovian 
memory effects. A regular method for the calculation of non-Markovian corrections to the Drude 
conductivity is presented. A quantitative theory of the recently discovered anomalous low-field 
magnetoresistance is developed for the system of two-dimensional electrons scattered by hard disks of 
radius a, randomly distributed with concentration n. For small magnetic fields the magentoresistance 
is found to be parabolic and inversely proportional to the gas parameter, 5p xx /p ~ — (uj c t) 2 /na 2 . 
In some interval of magnetic fields the magnetoresistance is shown to be linear Sp xx /p ~ ~ijJ c t in 
a good agreement with the experiment and numerical simulations. Magnetoresistance saturates for 
lo c t ^s> na 2 , when the anomalous memory effects are totally destroyed by the magnetic field. We 
also discuss magnetotransport at very low fields and show that at such fields magnetoresistance is 
determined by the trajectories having a long Lyapunov region. 

PACS numbers: 05.60.+W, 73.40.-c, 73.43.Qt, 73.50.Jt 



I. INTORDUCTION 

The problem of magnetoresistance in metal and semiconductor structures has been intensively dis- 
cussed in literature during the past three decades. A large number of both theoretical and experimental 
papers on this subject was published. Most of these works were devoted to the case of the degenerate 
two dimensional electron gas where the electrons move in the plane perpendicular to the magnetic field 
and scatter on a random impurity potential. In this situation only the electrons with energy close to 
the Fermi energy participate in conductance and the usual approach based on the Boltzmann equation 
leads to vanishing magnetoresistance. In other words, the longitudinal resistance p xx of the sample 
does not depend on the magnetic field B. This implies that the explanation of the experimentally 
observed B dependence of the longitudinal resistance should be sought beyond the Boltzmann theory. 

The intense exploration of this area began from the work of Altshuler et al^L where the experi- 
mentally observed in 2D metals and semiconductor structures negative magnetoresistance (MR), i.e. 
decreasing p xx with increasing B, was explained by quantum interference effects. It was shown that 
the magnetic field destroys the negative weak localization correction to the conductivity, thus resulting 
in decreasing longitudinal resistance. Since the first publication on the subject a vast amount of work 
has been devoted to its further exploration (see for review Ref.0). 

Two years prior to Ref. Q] there appeared a publication^ where a classical mechanism of negative 
magnetoresistance was discussed. The mechanism was investigated on the example of a gas of non- 
interacting electrons scattering on hard disks (antidots). It was shown that with increasing magnetic 
field there is an increasing number of closed electron orbits which avoid scatterers and therefore are 
not diffusive (see also recent discussions^^ of this mechanism). Electrons occupying these orbits 
do not participate in diffusion. As a result, the longitudinal resistance turns out to be proportional 
to the factor 1 — P, where P = exp(— 2tt/uj c t) is the probability of the existance of the circular 
closed orbit, which avoids scatterers (here lu c is the cyclotron frequency, r is the scattering time). 
Another classical mechanism was presented in Ref. 0, where the MR due to non-markovian dynamics 
of electrons trapped in some region of space was discussed. 

Notwithstanding these developments, the role of classical effects in magnetotransport was underap- 
preciated for a long time. A new boost to the research in this direction was given by Ref. (J, where it was 
shown that if electrons move in a smooth disorder potential and in a sufficiently strong magnetic fields a 
phenomenon called "classical localization" occurs. This phenomenon leads to the exponential suppres- 



sion of the longitudinal resistance: most electrons are trapped in localized equipotential trajectories 
and do not participate in diffusion. This work was followed by a series of W orksji2*ii*i2ii&i2ii£ii£ii! 
discussing different aspects of classical magnetotransport in 2D systems. It was shownii that for 
lower magnetic fields near the onset of the classical localization the magnetoresistance is positive, i.e. 
the longitudinal resistance grows with increasing magnetic field. In Refs. mil the combination of 
smooth disorder and strong scatterers (antidots) was considered. It was shown that in this system 
under certain conditions there are several regimes of the behavior of magnetoresistance depending on 
the strength of the magnetic field: first the longitudinal resistance decreases with growing field, then 
it saturates an d the n begins to grow. 

In R.efs.l.^4yi ^fil7y 8y9l lirillll2ll.'114l magnetor('sistaiK (' w as studied in a situation whore the magnetic 
field is classically strong, that is where the parameter (3 = lu c t is large. Recently, the region of 
classically small magnetic fields <C 1 was investigated numerically^i 6 - for the case of electrons 
scattering on strong scatterers. It was showni^ that memory effects due to double scattering of 
an electron on the same disk lead to a negative parabolic magnetoresistance (in the Ref. 0, where 
these processes were not taken into account, exponentially small MR was predicted). The numerical 
simulations 16 discovered a low-field classical anomaly of the MR. The anomaly was attributed to 
the memory effects specific for backscattering events. The simulations were performed for the 2D 
Lorenz gas which is a system of 2D electrons scattering on hard disks randomly distributed in plane 
with average concentration n. Magnetotransport in this system is characterized by two dimensionless 
parameters: j3 = w c t, and the gas parameter (3q — a/l — 2na 2 . Here a is the disk radius, ui c is 
the cyclotron frequency, r = l/vp is the mean free time and I ~ \/2na is the mean free path. 
The anomaly was observed in the case /3 <C 1, /3 -C 1. Both the numerical simulations and the 
qualitative considerations 16 indicated that at zero temperature the MR can be expressed in terms of 
a dimensionless function f(z) via 

Spxx - -M ( 4-) , (i) 



(2) 



p \Po 

where p is the resistivity for B — 0. Numerical results^ 6 - suggest that f(z) ~ z as z — > 0, yielding 

Spxx | | 
~ - \U C T. 

P 



The latter expression is in a very good agreement with experimental measurements of negative linear 
MR in a random antidot arrays^ 8 It is anomalous in two senses. First, it has a non-analytic depen- 
dence on the magnetic field. Second, it does not vanish in the limit of vanishing /3o, which is normally 
regarded as the expansion parameter for the corrections to the Drude-Boltzmann picture. This in- 
triguing behavior calls for a rigorous analytical theory of the effect, which would establish Eq. (JTJ and 
enable one to derive the analytical expression for function /. 

In this paper we present a detailed theor y o f the anomaly and give an expression for f(z) (the brief 
description of our results was given in Ref.UJ). We find that at some interval, 0.05 < z < 2, function 
f(z) is linear in agreement with numerical simulations and experiment, but at z — > crosses over to 
a quadratic dependence. Thus, for (3 — > 0, Eq. J3> yields Sp xx /p ~ — /3 2 /A>- The limit (3q — > should 
be taken with care. While the small (3 expansion seems to be singular as a function of (3q, the region 
of (3 where this expansion is valid shrinks as /3q — > 0. For z — > oo, / saturates at some constant value. 
Therefore, the full variation of Sp xx /p is of the order (3q. In other words, the anomalous MR is strong 
but it exists in a small region of magnetic fields. 



II. QUALITATIVE DISCUSSION OF THE PROBLEM 

The mechanism proposed in Ref. is linked to the memory effects arising in backscattering 
events. It has a close relation to the well known non-analyticity of the virial expansion of trans- 
port coefficientsji2i22i2ii22i2£ which we briefly recall. For B — the leading nonanalytic correction 
to resistivity, dp, is due to the processes of return to a scatterer after a single collision on another 
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FIG. 1: Backscattering process responsible for leading nonanalytic contribution to the resistivity at B — 
(a). For 6/0, the overlap area, Sb, between two corridors is small at large B (b). For <!> = 0, Sb decreases 
with B (c). For $ ^ and small B the values of Sb — So for time reversed trajectories have opposite signs 
(d,e). 



scatterer (see Fig. |^). The relative correction, Sp/p, is proportional to the corresponding backscat- 
tering probability, given by the product of e~ R / l d<&dR/l (which is the probability to reach scatterer 
2 without collision and scatter in the angle d<&) and the probability p to return without collisions 
from 2 to 1 (here I is the mean free pass). Assuming p — exp(— R/l) and integrating over intervals 
< <f> < a/R, a < R < oo, one obtain a 1 ^ 2 ?! 21 ! 22 ! 2 ? 

5p/p~ Jl \ d$e- 2 ^~/3 ln(l/2/3 ). (3) 
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In Ref. it was shown that the probability p is actually larger than exp(— R/l). Indeed, the 
exponent exp(— R/l) can be written as exp(— nS), where S — 2aR. It represents the probability of 
the existence of an empty corridor (free of the centers of the disk) of width 2a around the electron 
trajectory from 2 to 1. However, the passage of a particle from 1 to 2 ensures the absence of the 
disks centers in the region of width 2a around this part of trajectory (from 1 to 2) This reduces the 
scattering probability on the way back. The correct value of p can be estimated as 

P (R, $) = exp(-n(5 - So)) = exp(-i?/Z + nS (R, $)), (4) 

where 

S (R,$) = 2ai?-i? 2 |$|/2 (5) 

is the area of the overlap of the two corridors (see Fig.^i). For example, for $ = 0, we have So = 2aR 
and p = 1, which reflects the obvious fact that the particle cannot scatter, if it travels back along the 
same path. Taking into account the effect of "empty corridor" , we get 



r /R d$e-^ +nS °~p Q \n^y (6) 



6 JL ^ [°° — f a/R ^ ~-(2R/l)+nS _ a !„ ( C 
P Ja I 

where C is a constant of the order of unity. Thus, for B = the "empty corridor" effect simply 
changes the constant in the argument of the logarithm. 

The key idea suggested in Ref. was that for B ^ the area of the overlap of the two corridors, 
Sb, sharply depends on B, resulting in the observed MR. Indeed, it is seen from Fig. ^ that for 
P ^ Po, Sb — > 0, resulting in sharp negative MR 

6 Pxx r m r R ^ e „ 2R/l ^ nSB _ e „ So ^ (7) 
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The following qualitative explanation of the observed linear MR was presented in Ref. ^| The value 
ti(Sb — So) was estimated for <p = (see Fig. QJ;) to the first order in B as —nR 3 /R c = —R 3 /2alR c , 
where R c is the cyclotron radius. Assuming that this estimate also works at <j) ^ and expanding 
e ns B _ e ns to ^ ne fi rs ^ orc j er m q one g e t s Sp xx /p ~ —l/R c = —u c t. 

In fact, the physical picture of the phenomenon is more subtle. The contribution of any trajectory 
with <fi 7^ is cancelled to the first order in B by the contribution of the time-reversed trajectory, 
since the values of Sb — So are opposite for these paths (see Fig. 1 d, e). The cancellation does not 
occur only at very small <f> ~ 0. The integration in Eq. Q over 4> < (3 yields Sp xx / p ~ —(3 2 /(3 . Larger 
values of (f> also give a quadratic in (3 contribution to the MR. This contribution is positive and comes 
from the second order term in the expansion of e nSn - e nS ° in B. It follows from our results (see 
Eqs. I0,jn2j) that the contribution of small angles is dominant resulting in a negative parabolic MR. 
We find that the parabolic MR crosses over to linear at very small (3 w 0.05/3o, which explains why 
the parabolic MR was not seen in numerical simulations^ and experiment. 18 Note that for very small 
/3 < Pq contribution of the trajectories with a long Lyaponov region (Lyapunov trajectories) becomes 
important. In the Section ITTT1 we focus on the region (3 > /3q where the contribution of the Lyapunov 
trajectories is parametrically small. The Lyapunov trajectories will be discussed in Section ITVl 



III. CALCULATIONS 



A. Kinetic equation 

In this subsection we introduce the kinetic equation which is the starting point for the calculation 
of the diffusion coefficient and the resistivity of the Lorenz gas. 

We consider the Lorenz gas at zero temperature, assuming that the electrons participating in the 
conduction have the Fermi velocity vr. The diffusion coefficient D is given by 

D = IJ Q ^( v (0)v(t)> =\j dtj drdv(G fl )vv . (8) 

Here Gr = Gr(v, v , r, r , t) is the retarded Green's function of the Liouville equation and (...) 
stands for the averaging over the positions of the disks. The equation for Gr reads 

| + - ^ c [v x A] - Tj G R = S(r - r )*(v - v )J(t), 

G R = 0, for t < 0, (9) 

where vd/dr — lu c [v x d/dv] is the Liouville operator of the free motion in the magnetic field and T 
describes scattering on the disks. From Eq. © it follows that we need the time integral of Gr rather 
than the whole time-dependent function. This integral can be written as 
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/ dtG R (v,v ,r,r ,t) 
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S(v — vp) 
v F 



G Q , Q0 (r,r ), (10) 



where a, ao are the angles of velocities v and Vo respectively. Here we used the energy conservation, 
which implies that in scattering processes the absolute value of the velocity does not change. Using 
Eq. UlUfl . the diffusion coefficient Eq. @ can be rewritten as 

Vr f 

D = / drdrodada (G a!ao (r,r )) cos(a - a ). (11) 

Since (G Q) a (r, ro)) depends on a - ao and r — ro only, it is convenient to average over the position 
of the initial point and over the initial angle. Here S is the area of the sample. Integrating Eq. 
with respect to time the equation for G — G Q . ao (r, r ) is written as 



(e + L - f_ - f+) G = 5(r — r )S(a - ao), 



(12) 



V 



FIG. 2: Scattering of a particle on a hard disk. 



where 

L -n— - — — 

° <9r i? c da ' 

n = cos ae x + sinae^ is the unit vector in the direction of v and e -* 0. The interaction with disks 
is written in Eq. i|12|) in the form of collision integrali22i2£ The scattering operators T ± transform 
arbitrary function f(r,a) as follows, 

f +/(r, a) = J da'a{a - a')n+f(r, a'), 

f-f{r,a) = -f(r,a) / da' 'a(a - a')rT , (13) 



where 

± =n ± (r,a) = ^ <5(r - R, ± a). (14) 
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Here R.; are the positions of the disk's centers and 
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a((p) = ^ sin^ (15) 



2 



is the differential cross-section of one disk. The vector a depends on angles a, a' in the integrals 
Eq. H3J 

(n'-n)a 

a = a a . a > = = =, (16) 

\/2(l - n'n) 

and is pointing from the center of a disk to the scattering point at the disk surface (see Fig. |2J). 
Physically, operator T + describes influx of particles to velocity v at the point R^ — a, while operator 
T~ describes the outflux from velocity v at the point R^ + a. The averaged value of the n ± is equal 
to the disks concentartion 

(n+> = (n") = n. (17) 

The Bolzmann equation is obtained by averaging the Liouville equation with respect to the position 
of the scatterers and neglecting correlations. Indeed, using Eq. (|17|l one finds that the average collision 
operators Eq. l(T3|l are given by 



(T+)f(r,a)=n J da'a(a - a')f(r,a'), 

(f-) = -n [ da' a (a- a 1 ) = -l/l. (18) 



Replacing the collision operators in the Liouville equation Eq. (|12l) with their averages Eq. (|18fl one 
obtains the Boltzmann equation 

L Gl ao (r)-n J da'cr(a - a')(G B <ao (r) - G£ ao (r)) = 6(r)S(a - a ). (19) 

The Bolzmann-Drude diffusion coefficient Dq in the absence of the magnetic field is found from Eqs. 
(|f 9|) and (|f f ll as follows. Integrating (|19|l with respect to r and expanding the Green's function G B 
in angular harmonics one finds 

f Gl ao (v)dv = i-£ l m e im ^\ (20) 

m 

where 

^=n/^)(l-cos(m<,)), (21) 

In particular, 

ii = kr = y = «r- ( 22 ) 

is the transport length. Substituting Eqs. (|271| and (|2"2l in Eq. fTTJ) one finds 

A) = 4^ (23) 



B. Perturbative expansion 



In this subsection we derive the perturbative expansion for the average Green's function in order to 
take into account non-Markovian corrections which are absent in the Bolzmann-Drude picture. 
Introduce the operators 

Sf ± =f ± - (f ± ) (24) 

which describe the fluctuation of the collision integral with respect to its average value. Using these 
operators the formal solution of l|12|) . G = (e + Lq — T~ — T + ) _1 can be written as the following series 



(G)=G B + G B (5f»G B 6f u )G B + ••• , (25) 

where G B is defined by Eq. 1)19(1 . Here we took into account that {ST) = 0. The first term on the 
right hand side of Eq. (f^3|) gives the Drude-Boltzmann result described in the previous subsection. 
The rest of the terms in Eq. l)25ll provide a regular way for the calculation of correlations, which are 
absent in the Boltzmann picture. 

In the subsequent analysis of the perturbative expansion Eq. Il25[i we will extensively use the rep- 
resentation of the Boltzmann propagator Gb in terms of the ballistic propagator 

G" = (L + I/O -1 - (26) 

This is achieved by expanding the Boltzmann propagator G B = (Lq + 1/1— (T + )) _1 as a sum over 
the number of scattering events represented by the operator (T + ) 



G B = G~ + G~ (T+)G- + 



(27) 
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FIG. 3: Expansion of Boltzmann propagator G B in a series in ballistic propagators G 
N=l AT = 2, 3,... 
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FIG. 4: Diagrams, corresponding to the process shown in Fig.^i. Diagram (a) does not take into account the 
"empty corridor" effect and should be renormalized by diagrams (b). 



The diagrammatic representation of the expansion (|27() is shown in Fig. Q. Each cross in this figure 
corresponds to (T + ) and each solid line to a ballistic propagator G~ . For future reference we write 
the explicit expression for the ballistic Green's function in zero magnetic field and in weak magnetic 
fields. For B = the ballistic propagator conserves the velocity. In this case its kernel is given by 

G- Q0 (r) = S(a - a )S(a - a[r]) 6XP ( ~ r/0 . (28) 

Here a[r] is the angle of the vector r. Magnetic field rotates the velocity vector with the cyclotron 
frequency. For small magnetic fields, <C 1, we have 

G^ ao (r) « 6(a - a[r] - r/2R c )S(a - a[r] + r/2R c ) 6XP ( ~ r/0 . (29) 



C. Ballistic returns in perturbative expansion 



In this subsection we show how the processes of ballistic returns discussed in Section [n] arise in 
the perturbative expansion Eq. I|27|l . Consider the second term in Eq. (|25|l . This term describes the 
memory effect due to diffusive returns. As discussed in Section [H] the main contribution comes from 
returns after a single scattering. This process is described by the diagram shown in Fig. Up,. The 
dashed line corresponds to the pairings (ST^ST") (fj,,u — ±), external wavy lines to the diffusion 
propagators G B . The internal line corresponds to the Boltzmann propagator Eq. I|27|l truncated at 
one scattering G~ (T + )G~ . 

Four combinations of (±,±) at the ends of external dashed lines in the diagrams shown in FigQJi 
represent four different types of correlation at a given point R. To see this consider the pairing 
(ST^ST"). By virtue of Eqs. and 114|) this pairing is proportional to the density-density correlation 




FIG. 5: Electron trajectory related to the pairing (+,+). An electron scatters twice (at points ri and r2) on 
the surface of a disk, placed at point R. For a fixed position of the backscattering point (r, = 0), the process 
is parameterized by R, and by the scattering angles <fio,tpf 
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FIG. 6: Electron trajectory for the ( — , — ) pairing. Electron does not experience any real scattering except 
as scattering at r* = 0. The diagram ( — , — ) gives the probability to pass twice the region around point R 
without scattering. 




FIG. 7: "Unphysical" electron trajectory related to the pairing (—,+): electron scatters on a disk and after 
a while passes the region occupied by this disk without scattering. 

function ((^(ri, ai) — ra) (n^ (r 2 , a.2) ~~ n )) f° r /i = = ±- Assuming that disks are randomly 
distributed over the sample we get for these functions 

(K(n,ai)-n)K(r 2 ,a 2 )-n)) = (30) 

n J dR<5(r 2 — R + /1 ai)(5(ri — R + v a 2 ) 

It is natural to interpret the vector R as the position of the center of the disk on which a double 
scattering occurs. However, this interpretation is valid for diagram (+,+) only. As shown in Fig. 
this diagram corresponds to the situation where an electron experiences two real scattering processes 
on the disk placed at point R. The physical interpretation of other diagrams is more subtle. The 
diagram (— , — ) (see Fig. 0) does not correspond to any real scattering at point R. It just allows one 
to calculate correctly the probability for an electron to pass twice the region of the size a around 
point R without scattering. To interpret the diagrams (+, — ) and (— , +), note that in the Boltzmann 




FIG. 8: "Unphysical" electron trajectory related to the pairing (+,—): electron passes the region occupied by 
a disk without scattering and later on scatters on this disk. The diagram related to such trajectory corrects 
the Boltzmann result by subtracting the probability of this unphysical event. It is convinient to parameterize 
this trajectory by the scattering angles for the (+, +) process. 
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FIG. 9: Graphical representation of Eqs. 1311 . 132H 



picture, which neglects correlations, the following process is allowed. An electron scatters on a disk and 
later on passes through the region occupied by this disk without a scattering (see Fig. [7J (analogous 
consideration is valid for diagram shown in Fig. (SJ. The diagrams (+, — ) and (— ,+) correct the 
Boltzmann result by subtracting the contribution of such unphysical processes. 

As follows from the qualitative discussion Section [HI the processes of ballistic returns after a single 
scattering give rise to the leading (non-analytic) correction to the Drude-Boltzmann result in zero 
magnetic field. However, taking into account diagrams shown in Fig. 0Ji (describing processes shown 
in Figs. l5l8|l is not sufficient for calculation of the low field anomaly of the MR. Actually, diagrams in 
Fig. QJi do not contain the "empty corridor" effect. We will show that the correct description of the 
memory effects, specific for ballistic returns, requires the renormalization of diagrams^ by diagrams 
in which the internal dashed lines contains (— , — ) pairings only. Physically, the iV'th order diagram 
of the type^Js represents the N — l'th order term in the Taylor expansion of the exp(— uSq) in the 
qualitative estimate © . Such renormalization play the key role in the quantitative description of the 
anomalous MR. 

In the next subsection we derive analytical expression for diagrams QJi. Then we generalize the 
calculations to account for the corridor effect and calculation of anomalous MR. 



D. Non-analytical corrections to the zero-field resistance neglecting the "empty corridor" 

effect 

In this subsection we use the perturbative expansion Eq. (|25[1 to derive analytically the leading 
correction to the diffusion coefficient. This correction is due to the processes described by diagram 
Fig. The operator expression of the correction to the Green's function related to this diagram is 
given by 

G B 5T^G-{T+)G-8T V G B , 

where the underbracket stands for the pairings of operators 5T 11 and 8T U . Substituting this correction 
into Eq. iJHJ, and using Eqs. (|2L)|) . I|21I) to integrate over r, ro, a, ao we get 

SD = ^ / d0od0M 0A cos(/3 - ft) (31) 

where 

Mp,p = ^ Y] / dr^dr 1 dr 2 daida2d/3idP2 (32) 
S v*>=± 3 

(ST^(r 2 )G- 2t02 (r 2 - r*)na(ft - ft) G ft , Qi (r, - v 1 )6T^ o (r 1 )) . 



Here we used Eq. (|18f) for (T + ). Graphically, these calculations are presented in Fig. |5J The pairings 
(STg a2 (^2)5T^, i o (ri)), entering Eq. (|32|l . are given in Appendix A for four possible combinations of 
/i, v. The pairings only depend on the difference i"i — r 2 . Therefore in Eq. i|32|) one can remove the 
integral over r* and the sample area S in the denominator and put r* = in the integrand. In other 
words we put the origin of the coordinate system at the backscattering point r». Next we integrate 
Eq. IML'I) over the angles a±, a 2 , 0i, @2- This integration can be easily done taking into account angle 



dependent delta functions entering in G a2 g 2 {t 2 ) and Go ai (—Ti) (see Eq. ljz?8)) ). As a result we get 
ai = 0i = a[—ri], a 2 — 02 — ot[r 2 ] and 




In derivation of Eq. (|33[) we took into account that the pairing of collision operators {8T^ {yi)8T v {t 2 )) 
vanishes for |ri — ra| > 2a, which is clear from its physical meaning and also can be seen from the 
delta functions in the factor J(a, b) (see Appendix A). At the same time both r% and r 2 are of the 
order of the mean free path I a. Therefore, 02 — Pi — ^[^2] — ce[— ri] ~ tt, and a {02 — 0x) ~ 7r. 
Physically this means that the typical scattering angle in the processes of ballistic return are close to 
7T. Next, we use Eq. (|33H to evaluate the contributions of different pairings ji, v. 



(+>+) pairing 

Consider first the process in which \i = +, v = +. Substitute equations i|A3|l and l|33|) in Eq. (|31|l 
and integrate over , r 2 . It is convenient for our purposes to perform this integration for a fixed value 
of the vector R. Then the delta functions in Eq. (|A7|I ensure that ri and r 2 lie on the circle of the 
radius a centered in the point R 

n =R-a a[ _ ri ] )( g (34) 

r 2 = R - a.p,a[r 2 ]- 

The graphical solution of these equations is shown in Fig. 03 ( vectors na, ng in this figure are the unit 
vectors in the directions of and 0o respectively). Since a is small (o<() one has approximately 

ri*r 2 « R. (35) 

This accuracy is sufficient for the calculation of the diagram in Fig. 0Jl (since we calculate it in the 
lowest order in a/ 1 ). Upon integrating out the vectors ri, r 2 the integrand in Eq. (|3*l"1) depends on 
R. Finally we get 

6D ++ = v F lln 2 CT{ir) r dRdf3d(3o cos(/? _ (36) 
4tt J 

-2R/1 

a{0 - a[K])a{a[-K] - 0o)^~ 

Introducing new variables tpo = 0q — a[— R] and <ff = a[R] — and taking into account that a[R] = 
a[— R] + 7r and a{—tp) — cr{(p), the relative correction to the diffusion coefficient is written as 

SD^ ^ Tlltr f dFt —2R/] / \ / \ t \ 

dipodipf e < a{tpf)a{(po)cos((p + ipf) (37) 



Do 41 J R 

^From Fig. [S] it is seen that ifiQ, iff are the scattering angles. 



(+,-) pairing 

Next we present calculations for /1 = +, v = — . Like the calculation of the diagram (+,+), wc 
substitute Eq. I|33|l and Eq. i|A5() into Eq. (|31|l and integrate over ri, r 2 with the use of Eq. (|A7I) . 
The values of ri and r 2 should be found from the following equations 



r a — R + a Q [_ ri ] Q / 

r 2 = R a/3, Q [r 2 ]- 



(38) 



The graphical solution of these equation is presented in Fig. |S| After integrating out vectors ri , r 2 
and the angle 0q, we get 



SD +- = VFljn 2 Cj{TT) 
47T 

e -2R/l 

o{p - a[K])a(a[-K] - «')-^- 



J dUdfida' cos(/3 - a[-R]) (39) 



In derivation of this equation we have used Eq. I|35|) . It is convenient to rewrite this equation using 
as integration variables ipo, iff, which are the scattering angles for the process (+,+). They are 
expressed in terms of a' and j3 as (po = a[— R] — a' , tp/ = a[R] — (3. In these variables Eq. I|39|l reads 

SD+- rc/tr f dR 



= if / ¥^ ^ e-^W/Mw) cos(^) (40) 



It is seen from Eq. I|4U[1 (see also Fig. [SJ that the contribution of the process (+, — ) can be param- 
eterized by the angles <po and </?/, which are the scattering angles for the process (+, +). Analogous 
calculations can be easily done for other types of correlations. 

Summing the different contributions SD — 8D ++ + SD" 1 + SD ^ + 8D~~ we get the following 
expression for the diagram Fig. 0^, 

bp SD nl tl . f°° dR 2J?/; [ 21T , ( 2lx , , . , , 

tr -Re / —e- 2R/l / d<A) / dtp f <r(<p Q )a{<p f ) (41) 



P D 41 J a R 

(l-e^) ( l-e-) = f 

This equation is the exact expression for the non-analytic correctioni2i22i2ii22i2ii to the Drude- 
Boltzmann resistivity which was qualitatively given in Eq. |3fl. Four terms in the product (1 — 
e l w)(l — e l(ff ) = 1 — e lVo — e l(pf + e^ Vo+¥ ' / - ) correspond to four types of correlations discussed above. 



E. "Empty corridor" effects on the zero- field resistance 

In this subsection we use the perturbation theory for the quantitative derivation of the effect of the 
"empty corridor" discussed in section ITT1 

Ea. H41|) takes into account one pairing of operators 8T^ and 8T V . The terms containing N pairings 
(N dashed lines) are typically small as @q . However, there is a series of diagrams, shown in Fig.^J), 
whose contribution is of the order /3o (see Ref. l2lTl . The internal dashed lines in this series only 
contain pairings (5T~8T~). Below we show that the series Fig. 0Jd accounts for the effect of the 
"empty corridor". More precisely, we prove that the TV— th order term in this series corresponds to 
N — 1 term in the Taylor expansion of the exp(nS'o) in Eq. Q. 

The addition of the diagram Fig. 2t> to the diagram Fig. 0Jl leads to the following replacement in 
Eq. 

G o 2 , /32 (r 2 )G ft!tti (-r 1 ) - K%£( ri ,T 2 ), (42) 

where the function K^'^(rx,T2) can be found from the Dyson equation (See Fig. I1(J|I. Since the 
operator ST~ does not change the velocity angle, one can search the solution of the Dyson equation 
in the form 

K%%(t 1 ,t 2 )= (43) 
G a 2 , /32 (r 2 )G ft!Qi (-r 1 )r(r 1 ,r 2 ,<i>), 



r 2 , «2 1"2, «2 1"2, a-2 




FIG. 10: Dyson equation, describing the renormalization of backscattering events by the "empty corridor" 
effect. N— th order term in the these stair-like diagrams corresponds to the (N — 1)— th term in the Taylor 
expansion of exp(nSo) in Eq. @. 




FIG. 11: Pairings (ST~ST~) do not change the electron trajectory. Point r (see Fig. 1101 lies on the line 
segment connecting point ri and the origin and point r' lies on the line segment connecting point r-2 and the 
origin 



where $ is the difference between the backscattering angle and ir (See Fig. Hill) ). While substituting 
Eq. I|43|) in the Dyson equation (which is shown graphically in Fig. (|10f) . integrating over angles a, a', 
and using the identity 



5(a[r 2 - r] - a[r]) = ^ ^ 5{a[r] - a[r 2 ])6»(r 2 - r) (44) 

(and analogous identity for vectors r',ri )we get 

r(ri,r 2> $) = l+ (45) 
/ / drdv'^^^^^t^ (st-^ (r)6T- { _ Ti] (r')) T(r, r' , 

J r<r2 J r'<ri 11 

where ST~ (r) are defined in the Appendix Further calculations can be done in two different ways. 
It is useful from methodological point of view to discuss both of them. One way is to first integrate 
Eq. (|46|l over the angles of the vectors r and r' (see Fig. HUfl . The delta- functions entering in Eq. (|46|l 
ensure that point r (see Figs. HUIllfl lies on the line segment connecting point ri and the origin and 
point r' lies on the line segment connecting point r 2 and the origin. This allows to reduce the Dyson 
equation to the closed relation for function T(ri, r 2 , $) 

r(n,ra,$) = 1+ (46) 
I " dT ' £ ^ ( 6T ^ ] {r)ST ^ (r ' } ) r(r ' r '' $) 

where the pairing 

(5T- [r2] (r)5T- hri] (rO) « h(r, $)5(r - r'), (47) 

is calculated in Appendix 151 Here 

h(r, $) = (2a - |$| r) 9 (2a - |$| r) (48) 
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FIG. 12: Width of the overlap region as a function of a distance from the origin 



is shown in Figll2l As a result we have 

r(n,r 2 ,$) = 1+ (49) 

^min[ri ,r 2 ] 

drh(r,<&)T(r, r, $) 
This equation has an evident solution 

/ rmin[r 1 ,r 2 ] \ 

^(^l,^2,'I , ) = exp n / drh(r 7 §) . (50) 



Eq. 1)50(1 may be derived in an alternative way. This way allows to understand on the formal basis 
why one should renormalize the diagrams in Fig. 0Jl by the (5T~ 6T~) pairings only. Let us integrate 
in Eq. over r, r', using Eq. (|A4|) . Doing this, we get 

r(ri,r 2 ,$) = 1+ 

da'da" [ ^ R £(a[R + a a [ r2 ], a >] - a[r 2 ])6(a[-R - a a [_ ri ], a >/] - a[-ri]) 

J i?.<min[ri ,r 2 ] ^ 

a(a[r 2 ] - a>(a[ri] - a")T(R, R, $). (51) 

Here we took into account that r(r, r', $) is a slowly changing function of r, r' (it changes on a scale 
of the order of I). Since Eq. SKty provide |r - R| ~ |r' - R| ~ a, we put T(r, r' , $) w T(R, R, $). Next 
we make use of the identity 

J da'a(a ~ a')S(a[R + a«, a /] - a) w R 6(a/R -\a-a[R]\) (52) 

to integrate over a', a". After this integration we get 
r(n,r 2 ,$) = 1+ 

/•minfri ,7*2] 

n/ RdRda[R]6{a/R- \a[r 2 ] - a[R]\)9(a/R- |a[r x ] - a[R] |)r(i?, i?', $) (53) 

Jo 

Integration over a[R] leads again to Eq. 1)50(1 . 

^From Eq. ((35(1 it follows that we need to know T for n ~ ^2 ~ -R. Using Eq. 1(50(1 we find 



r(i?,i?, $) = exp drh(r,<S>) 



e nS (iJ,$) (54) 



where Sq(R, $) is the overlap between two corridors, given by Eq. JSJ. The n-th order term in the 
Taylor expansion of the exponential ((54(1 corresponds to a diagram in Fig. ^3 with n+ 1 dashed lines. 
Indeed, the n-th order term in the Taylor expansion contains n integration over coordinates of n 
scatterings of the type — , — . What remains to do to get the resistivity correction is to express $ in 
Eq. ((54|) via angles ipo ,(pf. To do this, the precision Eq. 1(35(1 is not sufficient. More specifically, we 



will need to know the angle between vectors ri and r2 to the order a/l. For this purpose it will be 
sufficient to make replacements a[— r\] sa a[— R] and a[v2\ ~ a[R] in the arguments of vectors a in 
Eq. Ij34|l . Then the angle $ is calculated as (see Figs. 



$ ~ ° [cos(^ /2) + cos(p//2)] (55) 
K 

This equation is valid for < tpo < 2n, < iff < 2n (since cos(</j/2) is not periodic with 2-7T one 
should specify the integration limits) . It worth noting that Eq. (133H ensures that the argument of the 
function in Eq. (|48(l is positive and the overlap width and overlap area are given by 



/i(r,$) = 2a- |$|r 



S (R,<f>) = h(r,$)dr = 2aR- |$| i? 2 / 2 
Jo 



(56) 



Summing the diagrams Fig. 0Jd together with Fig. one gets an exact equation 

roo dR ,2* z-Stt 



Sp 


nl ilT) 
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- e lV0 ){l - 



^-e 2R ' 1 I dipo I dtp f a(tp )a(tp f ) 



(i 



2/?o ( C_ 
2A) 



(57) 



instead of qualitative estimate Eq. @ . 

Here C w 1.8 is a numerical coefficient. Thus, addition of the series Fig. HJ) to Fig.0^, leads to the 
following renormalization: ln(l/2/3o) — ► ln(C/2/3o). 



F. Estimate for neglected diagrams 



In this subsection we use the derivation presented in the previous subsection to show how to select 
relevant diagrams. The derivation was based on identity l|52|) . The left hand side of this equation is 
proportional to the impurity cross section a(a — a') ~ a. However, the right hand side is parameteri- 
cally larger ~ R, provided that \a — a[R]\ < a/R. One can show that the function similar to the one 
on the right hand side l|52(l arises each time when one of the index fi or v in the pairing (ST^ST") 
is equal to (— ). (for a pairing (5T~5T~) we have two functions of the type l|52|l . see Eq. J53J)). For 
the case when fi or v equals to (+) we have a differential cross section er(a — a[R]) ~ a instead of 
function l|52l) . Now we are ready to estimate different types of diagrams. Insertion of one additional 
(— , — ) pairing into a diagram in Figs. 15151 gives two additional Green's functions and, consequently, 
multiplier ~ l/^ 2 , and two functions of the type (|52|) giving a multiplier ~ I 2 . One should also multiply 
on disk concentration n and integrate over dR. Due to 9 functions in Eq. I|52|) the integration area 
dR is of the order la. Combining all the multipliers together we have nla ~ 1. Therefore, addition of 
(— ,— ) pairing does not lead to any smallness. In contrast to this, insertion of (— ,+) pairing leads 
to smallness ~ a/l. Indeed, the only difference from the case of (— ,— ) pairing is that one should 
replace one of the function H52fl by the corresponding cross section. This leads to the change of one 
of the multiplier of the order I by multiplier of the order a. Following this line of reasoning one could 
conclude that insertion of (+,+) pairing leads to the smallness of the order of (a/l) 2 . Actually, the 
smallness arising from insertion of (+, +) pairing is of the order of a/l. Indeed, while the replacement 
of each of two function of the type 1521) by corresponding cross section leads to the relative smallness 
of the order (a/l) 2 , the integration area in this case is not restricted by angle dependent 6 functions. 
As a result, the integration area, dTL ~ I , is larger by a parameter l/a compared to the case of (— , — ) 
pairing. Note that estimates presented above do not work for the diagrams with one dashed line 
Figs. 15181 In this case one should integrate over initial and final scattering angles. One can easily see 
that integration of Eq. (|52|l over angle a leads to the multiplier of the order of a. This implies that 

in this case all four pairings, (++),(H — ),( — h), and( ) are of the same order, of a/l (see Eq. flU) )- 

Note that the higher order diagrams which are small in the parameter [3q may turn out to be relevant 
for the MR at very low magnetic field (see Section [lV|l . 



FIG. 13: The magnetic field changes the backscattering angle <f> = $ + R/R c . The dashed (solid) line in 
represents electron trajectory for B = (B / 0). The electron trajectory for B 7^ is parameterized by the 
scattering angles 5 = 



G. Anomalous magnetoresistance 



In this subsection we generalize our calculations for the B ^ case assuming that B is small 
((3 <C 1). The main contribution in this case still comes from diagrams in Figs. 0Ji,b. Consider for 
example diagram (H — h). Let us compare the process of double scattering described by this diagram for 
B^O (see Fig. ^J) with the same process for B = (see Fig.0). For fixed points ri, r 2 and r» = 0, 
one can see the following differences. First, the scattering angles ipo and cpj acquire small corrections 
of the order of R/R c ~ l/R c = ft *C 1. Second, the parts of the electron trajectory corresponding to 
free ends of the picture become curved. The backscattering angle (j> increases by the value R/R c 

cj> = $ + R/R c , (58) 

where $ is the value of backscattering angle for B = given by Eq. H55[l . Finally, the overlap area of 
the corridors changes because the trajectories become curved (see Fig. 1131) . The corrections to ipo and 
iff lead to small relative corrections to the resistivity of the order of flo(3 2 and can be neglected. The 
same reason allows one to neglect the curvature of the incoming and outgoing parts of the trajectory. 
Therefore, the only relevant difference is the change of the overlap area of the corridors. The solution 
of the Dyson equation is analogous to the B = case. The points of intermediate integration r and 
r' lie now at the segments of cyclotron circles (from ri to r» and from r* to r 2 ). The pairing of two 
operators ST~ is still given by Eq. with the replacement h(r) — > fts(r). For (3, /3q <C 1 the overlap 
width is calculated as 

h B {x) w (2a - \<frr - r 2 /R c \) 9 (2a - \<fyr - r 2 / ' R c \) , (59) 

where 9 is the Heaviside step function. Therefore, the only difference from Eq. I|57|l is that one should 
replace Sq — > Sb, where the overlap area is given by 

S B (R,(f>) = I drh B (r), (60) 



The value of Sp xx /p is obtained from Eq. I|57|l by replacing e ,lSo to e nS ' B — e nS ° 

(1 - e ilp0 ){\ - e lv f ) ( e »SB(R,*) _ e nS (R,t>)^ i ( 61 ) 



Introducing dimensionless variables T — R/l, z = fl/0o we get Eq. where function /(z) is given 
by 



cos 



Here 
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dip / d(p f 



(62) 
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cos(<^o/2) + cos(</?//2) zT 
C = 2y ^ "2"' s ° = Sz ^Q- 



Function /(z) has the following asymptotics 



0.32z 2 for z ^ 

0.39- 1.3/Vi for z -> 00. 

In the interval 0.05 < z < 2, /(z) can be well approximated by linear function 

f(z) w 0.032 (z - 0.04) for 0.05 < z < 2. 



(63) 



(64) 



(65) 



Next we discuss the parabolic asymptotics more carefully. We will show that in this asymptotic 
region there are two contributions of different signs to the magnetoresistance: a negative contribution 
coming from the trajectories with very small $ such that |$| < j3, and a positive contribution coming 
from larger angles. This considerations will be used in the next section, discussing trajectories with 
the long Lyapunov region. 

First, we write the difference of two exponents e nSB —e nSo as follows (e Tl ( SB-S °) — l)e TlSo and expand 
the equation in the bracket in the Taylor expansion up to the second order 



e nS B _ e nS ^ e nSo | nSS . 



n 2 6S 2 



(66) 



SS — Sb — Sq. 



We consider the case z <C 1 (/? <C /So)- O ne can easily see that in this case the expression 2a — 
\<jyr — r 2 /R c \ , entering in the argument of 8 function in the Eq. I|59|) is positive at all values of <f> 
(we take into account Eq. (|55|l and have in mind that r < I). Therefore, the difference SS(&) = 
Sb(&) — Sq(&) is expressed as 



<SS($) 



dr r(|$| - \<f>+(R-r)/R c \) 



6J?, C 3R C 



9(R + R c <f>) 



for $ > 
for $ < 0. 



(67) 



To sum the contributions of different electron trajectories, we take into account that the time reversed 
trajectories have the same statistical weight. Indeed, as seen from Figs. (|5I8|I . (|13|) . time reversion 
correspond to the change ipo — ► 27r — iff, ipf — > 27r — <pg. This transformation does not affect the 
factor a((fo)a(tpf)(l — e llpo )(l — e lv/ ), entering Eq. (|61|l . At the same time, the angle $ changes to 
— $ under this transformation (see Eq. (J5HJ). From Eq. I|67|) we find the variation of the overlap area 
averaged over two time reversed trajectories 



6S(9) + 8S(-$) 



n(R-R c m 

GRr 



6(R-R C \$\ 



(68) 



This expression is of the order of /3/A) for |$| < R/R c ~ B and is equal to zero for larger angles. 
We conclude therefore, that time reversed contributions do not cancel only in the region of small 
angles 3> < B. Since the total variation of $ is much larger, of the order of Bq (see Eq. we can 

replace the expression in the right hand side of Eq. Ij68|l by (5-function, writing n(8S($>) + SS(— 3>))/2 w 
— (/3 2 / 24/3 )(i? 4 /? 4 )(5( ( I ) ). Keeping the leading terms of the order of /3 2 only, we obtain 

n6S(<S>) , e n5S(-#) fl2 n4 «2 n6 



2 24/? Z 4 v ' (24/3 ) 2 Z 6 

The second term in this equation comes from the averaging of the quadratic term in the exponent 
expansion n 2 SS 2 /2 over two time reversed trajectories. This term is positive and partly compensates 
the negative contribution of small angles. Using Eqs. . . l|f)T ]l . after some algebra we get the 
low-field asymptotic of MR as the sum of the negative and positive contributions discussed above 

J^ = -f-+A^-^-0.32^-, 70 
P 2po Po Po 

where numerical coefficient A is given by 

„ f v C , ,n sin 2 a sin 2 8 cos(a + 8) „„„ 

A = 320 / / dadB—r, —^cp- w 0.18 (71 

Jo Jo 4 + cosa-cos/? 6 



IV. CONTRIBUTION OF TRAJECTORIES WITH LONG LYAPUNOV REGION TO 

MAGNETORESISTANCE 



The equations derived in the previous section give the contribution to the MR related to the 
processes shown in Figs. 115181) . the parabolic asymptotics lf7U|) starting to work when 8 becomes 
smaller than Bq. Such processes are related to the correlations specific for returns to the initial 
point after one scattering. As we show at this section for very low magnetic fields, 8 <C Bq, other 
correlations come into play. Specifically, we consider the contribution to the MR of the trajectories 
containing long Lyapunov regions. Such trajectories consist of the direct and the return paths and 
involve real double scatterings on some number of disks as shown in Fig. 1141 The divergency between 
the direct and the return paths is characterized by the Lyapunov length. Wc will call such trajectories 
"Lyapunov trajectories". In the diagrammatic series they are presented by the sum of the diagrams 
shown in Fig.^J Just as in the case discussed above, four different pairings are allowed at the ends of 
external dashed lines: (— , — ), (+, — ), (— , +) and (+, +). However, in contrast to the diagrams shown 
in Fig. 01 internal lines of "Lyapunov diagrams" contain pairings of (+, +) type as well as of (— , — ) 
type. Physically, (+, +) pairings corresponds to real double scatterings in the Lyapunov region. We 
will count such diagrams by the number N of correlated links in the Lyapunov region. The Lyapunov 
trajectory shown in Fig. 1141 corresponds to iV = 3. The diagrams in Figs. 01 discussed in Section ITTT1 
present a particular case of diagrams with Lyapunov region, corresponding to iV = 1. The correlated 
links of any Lyapunov trajectory are renormalized by the pairings (— , — ) as was discussed already for 
N = 1. The contribution of any diagram of such type to the resistivity is small as Bq . However, as will 
be shown below, these diagrams have a sharp dependence on the magnetic field at very small fields. 
We will show that for 8 < 8$ the diagram of the N-th order gives a contribution to the parabolic MR 
of the same order as the contribution of the diagrams with N — 1 already calculated above. 

Consider the Lyapunov trajectory with N correlated links. Denote by $0j 3>i, ■ ■ ■ , the angles 
between the segments of the direct and the return paths and by ipo, . . . , </3jv-i the scattering angles 
between successive correlated links as shown in Fig.^J The contribution of such a process has a sharp 
dependence on the magnetic field due to the magnetic field dependence of the overlap of the corridors 
surrounding the direct and the return path. This dependence is different for different segments of the 
trajectory due to the difference of the angles $„. Indeed as we have seen in the previous section if 
$ is the typical angle between the direct and the return paths then the characteristic scale for the 
magnetic field dependence of the overlap of the corridors is 8 ~ <&. As one can see from Eq. (|55f) . the 



typical value of the angle is of the order (3q. The smallest angle <&n ~ (5q corresponds to the last 
segment of the trajectory. Therefore it is this segment that should lead to the sharpest dependence of 
the resistivity on the magnetic field. Consider the contribution of the last segment at small magnetic 
field (3 <C Pq . For typical trajectories with <&jv ~ (3q there is a cancellation of the contributions of the 
time reversed paths to the MR (see discussion in the previous section). However, for a small fraction 
of the trajectories with $jy < /3 the contribution of the time reversed paths does not cancel and is 
proportional to the change of the overlap area n(55'($jv) + 5S(— <J?at))/2 ~ (3//3o- The phase space of 
such trajectories is proportional to (3/(5$ . Thus, the contribution of the Lyapunov trajectory with TV 
links to the resistivity in the region of small magnetic field is given by 



Po ffi ~ Pa 



for (3<f5$ (72) 



where the factor (5q is due to the iV-th order of the corresponding diagram. As we show below, the 
coefficient in this equation does not depend on the diagram order. It equals 1/2, thus coinciding with 
the coefficient in the small angles contribution to parabolic asymptotic for N = 1 (see Eq. Q7U))). It 
worth noting that only in the case N = 1, one should take into account the contribution of large 
angles, presented by the second term in Eq. 1)70(1 . For N > 1 such contribution is parametrically 
small. 

Before presenting rigorous derivation of above statements, we consider qualitatively the contribution 
of N-th order diagram for larger magnetic fields (3 > (3q . At such fields the first order contribution 
of the time reversed paths to the corridor effect does not cancel for trajectories with any values of 4> 
up to the maximal value $ ~ (3q . Thus, the factor (3/(5$ , which counts the phase space of relevant 
trajectories saturates at some constant value of the order of unity and one has 

r N 

^--/S/^- 1 for (3>/3» (73) 

Consider now the MR in the interval (5q +1 < (5 < (3q . In this interval the Lyapunov trajectories 
with the number of double scatterings smaller than TV give a parabolic contribution Eq. 1(72(1 to the 
MR. The trajectories with the number of double scatterings bigger than n give a linear contribution 
Eq. ((73(1 . the main contribution (5$ (3 coming from the trajectories with N + 1 double scattering. As 
a result one has 

^ <»-(%- A) for (5^<(5<(5». (74) 

P \ 2 J Pa 

Here CV ~ 1 and A is given by Eq. ((71(1 . For N 3> 1 the second term in the right hand side of Eq. ((74(1 
can be neglected and we get finally 



S Pxx ^ _N[P_ 
P 2 O ' 



for (5^ +i <P<^\ IV >1. (75) 



Equation Eq. 1(74(1 indicates that at very low magnetic fields, which correspond to large N ~ ln/3/ ln/3o 
the results of section IIII CI become incorrect and logarithmic renormalization of the parabolic MR 
occurs 

Sp xx ^ m/3 (5 2 
p ~ 21n/3 /? 1 ' 

Next we present a rigorous derivation of the contribution of the the 7V-th order Lyapunov trajectory 
for (3 -C (3q . As follows from qualitative considerations presented above, -C 1 for any < n < N. 
Indeed, <!>„ ~ (5q^o- Since the minimal angle, $at is on the order of (3, the maximal angle is small 
compared to unity, $o ~ (5/(5q 1. Therefore, the direct and return paths of the relevant Lyapunov 
trajectories are very close to each other. This allows us to characterize a scattering between n—1 and 
n links of the trajectory by one scattering angle ip n instead of two different scattering angles for direct 



FIG. 14: Electron trajectory with a Lyapunov region, containing N = 3 correlated links. 



and return paths. The dependence on the magnetic field comes from the last link of the trajectory. 
Expanding e n ' SB_s o ' — 1, like it was done in the previous section, we only keep the first term 

n SS{*s) + 8Si-* N ) a _jP_F&. (77 ) 

To write down the analytical expression for iV-th order contribution we should take into account for 
combinations of (±, ±) at the ends of diagram in Fig. (|15fl . This correspond to four types of scattering 
on the impurity in the Fig. (|14|) and leads to appearing of the factor (1 — e lipi )(1 — e llfl ) = |1 — e llpi \ 2 . 
Here we took into account that direct and return paths coincides in the first approximation and, as a 
consequence, scattering angles cpi, ip[ for direct and return paths are related to each other as follows 
ip± = 2tt — (p[ . The analytical expression for MR is given by 



— = -— J *pid*od Rt ■ -d Ks\l - e | ^ . . . ^ Woir 5(Z N ). 

(78) 

Here we integrate over distances between disks instead of integration over disks positions. The factors 
a 2 (ip n ) comes from double scatterings, the factor i2„ from denominators of the two Green's functions, 
describing the propagation along n— th link on the direct and return way. In writing Eq. I|78|) we also 
took into account the renormalization of correlated links by the (— , — ) pairings. As a consequence 
of this renormalization n— th link, which is passed twice (on the direct and return paths) comes with 
the factor e~ Rn / 1 instead of e - 2R ™/ 1 . What remains to do is to express final angle via initial angle 
<I>o- From simple geometric considerations we get 

^-^-..-^v. (79) 

Combining Eqs. (|78|l and 1)79(1 and performing the integral one gets 

5p N 1 B 2 

— f or N>1. (80) 

P 2 fa 



V. INTERPRETATION OF THE RESISTIVITY CORRECTION IN TERMS OF SMALL 
CHANGE OF THE EFFECTIVE SCATTERING CROSS-SECTION 

In our calculations of resistivity corrections we used Eq. JHJl as a starting point. Here we briefly 
discuss an alternative approach based on the accounting of the memory effects in terms of small change 



FIG. 15: The diagrams corresponding to electron trajectories with long Lyapunov regions 



of the effective cross-section. Consider first the (+, +) ballistic returns with given values of the angles 
ipo and ipf. The contribution of such processes to the resistivity (both for B = and for B ^ 0) 
contains the factor Re e^ vo+¥ ' / ^ = cos(<po + ¥>/)• Let us introduce now the scattering angle tp for the 
process (+, +) considered as a single scattering on a complex scatterer. For small values of ipo and ipf 
this angle is evidently given by(p = 7r + ipo + ipf. For arbitrary values of ipo, ipf the scattering angle 
reads 

f = V<po,<p s = (tt + <Po + ¥>/)(mod 2ir). (81) 

Such definition ensures that < tp < 2-7T. Next we introduce the integration over this angle: 
J dipodcpf cos(<^o + V 5 /) ■ • ■ = (—1) f dtpcosip J dipod<p fS(<p — <p Vo ,tp f ) ' ' ' ■ Analogous expressions can 
be easily written for the processes (— , — ), (+, — ), and (— , +). As a result, one can write the contribu- 
tion of four possible types of correlations as follows 



Re / dipo<kpf(l - e ltpa ){\ - e^f) ■ ■ ■ = 

dipcosip / dipodtpf (5(tp-ip Vo , Vf ) +8{tp - it) - 5(ip- ip vo , ) - 5(<p - <Po,<p f )) •■• • (82) 



The (— , — ) process is represented in Eq. I|82l) by 5(ip — it). Here we neglected a small difference (on the 
order of (3o) between the angle ip and it. By virtue of Eq. I|82|l one can easily check that the memory 
effects related to four types of ballistic returns can be accounted quantitatively in the frame of usual 
Boltzmann equation. One should just replace the scattering cross-section on one disk cr(<p) by the 
effective cross-section 

o-off(y) = <r(<p) + 5<tb(<p), (83) 
where S(tb(<p) is the field-dependent correction given by 

' - 2R l l I d,n„ I d l n, r r(,n^rr(,n^ e nSB 



Svb(v) =j { j -g-e 2H/t I dtpo I dip f (T(tpo)a(ip f ) 



Jo 



x [6(<p - <p<p 0tVf ) + 5{ip - 7r) - 8{ip - (p Vo ,o) - 5{ip - (po, Vf )}. (84) 
This correction does not change the total cross-section 

dip 5<7 B {<p) = 0. (85) 

In other words, the enhancement of cross-section caused by processes (+, +) and (— , — ) is accompanied 
by the reduction of the scattering due to the (+, — ) and (— , +) correlations. The resistivity correction 
is proportional to the change of the inverse transport length 

/ 1 \ f 2lT f 2 * 
8 [ j— 1=71 / d(p(l — cost/?) 5ctb(<p) = — n I d(pcostp 5<jb(<p)- (86) 
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FIG. 16: The value ol Sp xx /p/3o from Eqs. Q, 16211 (solid line) shown as a function of /3/Po together with 
the results of numerical simulations— presented for different values of /3o (triangles for /3o = 0.09 , boxes 
for /3o = 0.06 , circles for (3q = 0.03). Data for all numerical curves are shown for (5 < 0.3. Inset: The 
crossover from quadratic to a linear dependence at /3//3o ~ 0.05. This crossover was not resolved in numerical 
simulations. 



Using Eqs. (82), (JHUl and (JHSJ one can easily get the Eq. JHH for MR. 

Finally we note that this approach is easily generalized for calculation of the contribution of the 
trajectories having long Lyapunov region. For such trajectories one should replace the angles ipo and 
(pf in Eq. |H2by <pi and ip[ w 2ir — tpi (see Section Hvjl . For very small (3 such that (3q +1 < f3 < (3q 
and N 3> 1, the magnetic field induced correction to the scattering cross-section is expressed in a 
simple form 



6a B (<p) - 6a (tp) _ 2/3 2 ln/3 / _ cos 2 (y>/2) \ 
2a ~ 9/3 ln/3 V w / 



(87) 



VI. DISCUSSION 



In the previous sections we derived analytical theory of the low-field anomaly in the magnetore- 
sistance, caused by sharp dependence of the memory effects specific for backscattering events on the 
magnetic field. Next we compare our calculations with the results of simulations and with experiment. 
Note first that there is a parametrically small nonanomalous correction to Eq. due to returns after 
multiple scatterings. This correction is given by^ 

Sp' xx /p^-0.2(3 p 2 . 

To compare the results of simulations^ with the theoretical results in a wider region of parameters 
(3, /?o, we substract Sp' xx /p from the numerical curves. Theoretical and numericali& results are plotted 
in Fig. l(TB|) . in the universal units, Sp xx /pf3o versus z = /5//3 - It is seen, that the theoretical and 
numerical results are in a very good agreement. Experimental measurement of the MR in the system 
of antidot arrays^ agrees qualitatevely with our predictions. The experimentally observed MR was 
linear in magnetic field in a wide region of fields. This region corresponds to the interval of magnetic 
fields, where our results can be aproximated by linear dependence (see Eq. ©). The magnetic fields 
used in experiment were relatively strong, and the parabolic asymptotic was not achieved. The 
quantitative comparison with the experiment*^ is more difficult due to several reasons. First of all, 
in the structures, used in experiment, besides antidots there were also short range scatterers. Our 
preliminary estimates show that accounting for a short range disorder can change the results of the 



calculations. Secondly, the antidot distribution was not fully chaotic in the experiment. To be more 
specific, the antidots were randomly moved from the regular square lattice distribution by shifts on the 
order of 30-40 % of the lattice constant. Finally, the antidots did not have equal sizes, the uncertainty 
of the size being on the order of 50 %. In spite of this, a very good quantitative agreement with the 
experiment can be achieved by appropriate choice of the antidot size a in the uncertainty interval (see 
also discussion of the experiment in Ref. Note also that generalization of the theory for the case 
when antidot sizes vary is straightforward. What has to be done is to average pairings (8T^8T V ) with 
the distribution function for the antidots size s. This leads to the following modification of the obtained 
formulae: one should change a(ipo)a(ip /) — > a(ipo)cr(tp /) in Eq. (|57|l and hs — + h-B in Eq. I|59|l . Here 
over line means averaging over antidot sizes. We do not present here corresponding calculations, since 
the distribution function of antidot sizes is not known for experiment^ 

Next we briefly discuss several interesting unsolved problems. Note first that above consideration 
of Lyapunov trajectories is valid, provided that N < l//3o- Indeed, while calculating the N— th order 
diagrams we neglected small contributions of (+,— ) and (— ,+) pairings on the all correlated links. 
Accounting of such pairings is rather tricky and is out of scope of this paper. We expect, however, that 
such correlations may lead to a factor on the order of (1 — c(3q) n ( c ~ 1) in Eq. Therefore, for 

N > 1/Pq the contribution of Lyapunov diagrams might become N dependent. This gives a low-field 
limit for our theory 

(3X3^°. (88) 

We conclude that the behavior of the MR in the limit (3 — > remains so far unclear. 

In our calculations we fully neglected quantum effects. Such effects should decrease the "effect of 
empty corridor" due to the diffraction on the edges of the disks. In the situations where the magnetic 
field is not very small, (3 > (3q, the neglecting of the diffraction effects is justified if a ^> \/Xpl 
(Xf is a Fermi wavelength). This criterion ensures that diffraction effects on the edges of the disks 
are not relevant at the scales of the order of I. In the opposite case, a <C \fXpl, the diffraction 
should destroy the "corridor effect" , thus suppressing the anomalous MR. For small magnetic fields, 
when /3 <C /3q> the criterion for negligibility of the diffraction becomes stronger, a ^S> ^/XfL, where 
L ~ Zln(l//3)/ln(l//?o) is the characteristic size of the Lyapunov region. Another quantum effect 
which can be especially important from the point of view of the possible experimental realizations is 
the weak localization phenomena. The weak localization correction to the conductivity also has an 
anomalous dependence on the magnetic field. Moreover, the interpretation of the memory effects in 
terms of small change of effective cross-section discussed in Section^is very close to the interpretation 
of weak localization phenomenon developed in Ref.|2^. Similar to discussion in Section[3 the coherent 
enhancement of backscattering amplitude caused by weak localization, is accompanied by reduction 
of coherent scattering in other directions, the total cross-section do not change. The competition 
of the non-Markovian and the weak localization effects might result in new interesting phenomena. 
The study of such competition in a system with spin-orbit interaction can be especially interesting, 
because, in contrast to the weak localization correction, the "corridor effect" is not very sensitive to 
spin-orbit coupling. The detailed analysis of the quantum effects is a challenging problem which will 
be addressed elsewhere. 

We did not investigate the temperature dependence of the phenomenon. This dependence is related 
to the scattering by phonons (or electron-electron scattering) neglected in our calculations. It worth 
noting that importance of the electron-phonon scattering is expected to increase with decreasing 
magnetic field. Indeed, the potential of the electron-phonon interaction depends on time, therefore 
restricting the maximal length of a trajectory with the Lyapunov region L < L p h, where L p h ~ VFi~ p h 
and T p h is the temperature dependent characteristic time of the electron-phonon scattering. This 
implies that at small magnetic fields one should replace the logarithmic factor ln(/3)/ln(/?o), entering 
Eq. H7(jfl by a temperature dependent coefficient of the order of L p h/l- 



VII. SUMMARY 



We propose a theory of the negative anomalous MR in a system with strong scatterers. It is shown 
that the anomaly in the MR arises due to suppression of "empty corridor effect" by magnetic field. 
A detailed description of different types of non-Markovian correlations related to ballistic returns is 
presented. A method of diagrammatic expansion of the Liouville equation is developed which allows us 
to describe analytically the effects of "empty corridor" on ballistic returns. The analytical expressions 
for anomalous MR in different intervals of magnetic fields are derived [see Eqs. Q, Ij62|l . (|64(l . (|74Jl ]. 
The MR at very low magnetic fields was shown to be determined by the contribution of electron 
trajectories having long Lyapunov region. An interpretation of the memory effects in terms of small 
change of the effective scattering cross-section is discussed. The analytical results are shown to be in 
very good agreement with the numerical simulations and experiment. 
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APPENDIX A: 

This Appendix contains the explicit expressions for the pairings (ST^ST"). Using Eqs. (|13fl . (1 1 81) 
and (|24p. kernels of operators ST + , 5T~, entering in Eq. I1-S2H can be written as 

ST+ a , =a{a-a')(n + -n) (Al) 
ST- a , = S(a - a')ST- 

where 

ST- = ( da"a{a - a")(n~ - n) (A2) 



The functions n + and n in Eq. I|A1|I depend on vectors a(n',n) and a(n",n) correspondingly 
(see Eqs. (|14H . <|lfci|) ). The pairings of operators 8T^ and ST" can be calculated with the use of 
Eqs. l|5U |l .l|A"i" )l .l|A"2 |l . As a result we have 

<^+ Q2 (r 2 )^r+ i/3o (r 1 )) = a( / 3-a 2 )a(a 1 - / 3 )J(-a^ Q2 ,-a Ql ^ ) (A3) 
(«7(ra)*2£(ri)) = J da'da"a(p - a')a(0 o - a")J(^,*',*p a ,*") (A4) 
(ST+a^ST-in)) = tr{p - a 2 ) J da'a(fi Q - a')J{- a ^ a2 ,^ , al ) (A5) 



(STp{ ra )8T+ ufh (T 1 ))=<r(ai-0 o ) I da'a(0 - a') J( a/W , -a Qli/3o ). (A6) 
Here function J emerges as a result of pairings of — n)(n u — n)) (see Eq. (|3(J|) ) and is given by 

J(a, b) = n / dR£(r 2 R - a)<5(r x R b) (A7) 



APPENDIX B: 

To calculate /<5T^ r3 j(r)<5T^_ ri j(r')\ we first do the integral in the right hand side of Eq. (|A2|) . 

Integral of the term containing average density n is trivial and reduces to l/l. In the integral containing 
n~ it is convenient to introduce orthogonal coordinates x and y such that x axis is collinear with the 
velocity of the incident particle and write the delta function entering the definition (|14(l of n~ as 
S(r — Hi — a) = S(y — R? — a y )5(x — Rf — a x ). We proceed by lifting the y dependent delta function 
by integration over the angle a'. It is clear that the integral is vanishing for \y — R?\ > a. To calculate 
the integral for smaller a note that by virtue of Eq. (|16fl \da y /da'\ = |rf(a x n)/da'\ = o(ot — a') 
and therefore the integral of the y dependent delta function cancels the scattering cross-section in 
Eq. I|A2(1 . As a result we get 

6T-(r) =l/l-J2 S( - x ~ ~ a *) * Kl l l ~Y. 5{ - x ~ R i) 6i ( B1 ) 

i i 

Here Qi is a unit step function which is equal to unity when \y — R\\ < a and vanishes otherwise. 
Using Eq. IjBljl we can write for backscattering angles close to it (<!> ~ a/l) 

( ST a[r 2 ](r) ST a[- ri] (*')) « nh (%, *)8(x - x') » nh(r, $)5(r - /), (B2) 

where 



is shown in Fig. ^] 



h(r, $) = (2a - |$| r) 9 (2a - |$| r) 



